Formation of a Heavy-Fermion State in the 2D Periodic Anderson Model 
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We study the formation of a heavy-fermion state in the 2D periodic Anerson model. For U = 2, the 
density of states, imaginary part of the self-energy and effective magnetic moment all indicate the 
Kondo screening of local f electrons, leading to a coherent heavy-fermion state. For 17 = 3 and 4, 
the dominance of RKKY interaction over Kondo screening at low temperatures indicates a magnetic 
instability at zero temperature. A partial screening of magnetic moments, however, still gives rise 
to a relatively sharp peak at the Fermi energy in the density of states. 

PACS numbers; 71.10.Fd, 71.27.-fa, 72.15.Qm 
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Since the discovery of the heavy-fermion materials in 
rare-earth or actinide elements , a possible dramatic 
transformation of localized / electrons to a coherent 
heavy-fermion state has received considerable attention. 
Experiments in these materials have exhibited various 
unusual properties such as the huge coefficient in the 
magnetic susceptibility, specific heat and so on, which 
all indicate the formation of a coherent state with a large 
effective mass. The periodic Anderson model has 
been considered as the most promising candidate which 
might be able to describe the anomalous features of these 
materials. The Hamiltonian for the 2D periodic Ander- 
son model is 

{i,j)a iff 
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Here t is the hopping parameter for c electrons, V the 
hybridization energy between c and / electrons, and e/ 
the energy level for / electrons. U is the Coulomb repul- 
sion energy for / electrons, — {fi^fia-) and fi is the 
chemical potential controlling the total electron concen- 
tration. CiCT and fia- are annihilation operators for c and 
/ electrons at site i with spin a, respectively. When the 
hybridization term is eliminated by means of the equa- 
tions of motion for c and / electrons, the correspond- 
ing interacting Green's functions, Gc(k,uj) and Gf{k,uj), 
are expressed in terms of the self-energy for / electrons 



Gfik,L0) 



1 



UJ 



y2 



1 



y2 



(2) 



hopping. As long as / electrons are concerned, the 
interacting Green's function has the same structure as 
that for the usual Hubbard model with the inverse of 
the noninteracting propagator [G° (fc, w)]"-^ replaced by 

Lo^e f+u — : — . For later use, we also show the nonin- 

teracting {U = 0) energy dispersion due to hybridization, 
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where = — 2t(cos 

kx ~t~ cos ky) for nearest neighbor 



Recently the infinite and two j^] dimensional peri- 
odic Anderson models for the symmetric case have been 
studied by using quantum Monte Carlo (QMC) calcula- 
tions. In these works for a half-filled band, the authors 
found a strong competition between the Kondo effect 
and antiferromagnetic order depending on the size 
of V and U. Because of its symmetric nature, however, 
the model can not describe any possible formation of a 
coherent heavy-fermion state from localized / electrons. 
The 2D asymmetric periodic Anderson model has been 
investigated by McQueen et al. H] within ffuctuation ex- 
change (FLEX) approximation. These authors found a 
developing density of states near the Fermi energy with 
decreasing temperature, which was argued as an indica- 
tion of the Kondo resonance. However, they could not 
observe a Fermi liquid behavior in the calculations of the 
self-energy. This is possibly due to the insufficient treat- 
ment of strong local correlations and 2D spin fluctuations 
in the FLEX approximation [l0|] . 

Recently we have developed a new approximation 
scheme for the 2D Hubbard model which shows the es- 
sential features of the model This theory satisfies 
simultaneously the correct atomic limit for large frequen- 
cies (refiected as the correct asymptotic behavior of the 
self-energy at large oj) as well as 2D spin fluctuations 
due to the Mermin- Wagner theorem These two fea- 
tures are very important in the study of the 2D periodic 
Anderson model, in that the screening of local moments 
by conduction electrons is properly described by the first 
property and any possible magnetic instability at zero 
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temperature by the second. One more good feature of 
our formulation over FLEX approximation or second- 
order (in U) perturbation study is that the dynamical 
spin susceptibility can be computed accurately by im- 
posing the exact sumrules. In this Letter we present the 
strong numerical evidence for the formation of a coherent 
heavy-fermion state from localized / electrons by show- 
ing the density of states, low frequency behavior of the 
self-energy, effective magnetic moment and quasiparticle 
residue. 

In order to compute properly the spin, charge, and 
particle-particle susceptibilities which govern the inter- 
actions between electrons, we impose the following three 
exact sumrules to them |ll|Jlo| : 
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T and N are the absolute temperature and number of 
lattice sites, g is a compact notation for (q, ivn) where 
ivn are either Fermionic or Bosonic Matsubara frequen- 
cies. The dynamical spin, charge and particle-particle 
susceptibilities are calculated by 
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Xphil) ^^'^ Xppiq) irreducible particle-hole and 

particle-particle susceptibilities, respectively, which are 
computed from 
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where G'j{k) is the noninteracting Green's function for 
/ electrons defined earlier. C/jp, Uch^ and Upp in Eq. |^ 
are renormalized interaction constants for each chan- 
nel which are calculated self-consistcntly by making an 
ansatz [/.,„ = U {m^n,]) / {{ni^){n\)) Q in Eq. ^. By defin- 
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ing Uap, Uch, and Upp this way, the Mcrmin- Wagner the- 
orem as well as the correct atomic limit for large uj are 
satisfied simultaneously [0. Throughout the calcula- 
tions we fixed £f — 0.45, V — 1 and the total electron 
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FIG. 1. (a) Density of states and (b) imaginary part of 
the self-energy at the Fermi surface for U — 2, e/ = 0.45 
and V = 1 with decreasing temperature. The long-dashed, 
dashed, dotted and solid curves denote T — 1/16, 1/64, 1/256 
and 1/1024, respectively. The inset in (a) shows the density 
of states for T = 1/1024 in the wide range of frequency axis. 



concentration at 2.25. The unit of energy is t and all en- 
ergies are measured from the chemical potential fi. We 
used a 64 X 64 lattice in momentum space and performed 
the calculations by means of well-established fast Fourier 
transforms (FFT). It should be also noted that we used 
a real frequency formulation in Eqs. (|^)-(|^) to avoid any 
possible uncertainties associated with numerical analyti- 
cal continuation. 

We start in Fig. Q by studying the density of 
states (Fig. |l|(a)) and imaginary part of the self-energy 
(Fig. pb)) for U = 2 with decreasing temperature from 
T = 1/16 to 1/1024. As the temperature is decreased, a 
sharp peak develops at the Fermi energy in Fig. |^(a). At 
T = 1/16 this peak is completely absent. To find that 
this feature is associated with the onset of the Fermi liq- 
uid, we calculated the imaginary part of the self-energy 
in Fig. |l|(b). With decreasing temperature, the scattering 
rates progressively decrease and vanish at the Fermi en- 
ergy much faster than linearly in frequency. The log-log 
plot of the scattering rates vs. frequency for T = 1/1024 
shows ImYi{kp, uj) ^ w^'^^ in a relatively narrow interval 
of ±[0.005 - 0.025]. 

In order to establish a more firm ground that this is 
indeed due to the screening of magnetic moments by con- 
duction electrons, we show the effective magnetic mo- 
ment defined as Txsp{^, 0) upon decreasing the temper- 
ature (filled circles in Fig. |(a)). Below T = 1/128, the 
effective magnetic moment vanishes linearly in tempera- 
ture, a clear indication of the Kondo screening. Above 
this temperature, it deviates significantly from a straight 
line and appears gradually saturating at high tempera- 
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FIG. 2. (a) EfTective magnetic moment (filled circles) and 
quasiparticle residue (open circles), and (b) quasiparticle dis- 
persion (small filled circles) for T = 1/256 against nonin- 
teracting dispersion (dashed curve) in the upper branch for 
U = 2, e f = 0.45 and V = 1. The solid line in (a) is a linear 
interpolation at low temperatures. 



tures. The quasiparticle residue at the Fermi surface 
along the direction (0, 0) — (vr, tt) is also presented as the 
open circles in Fig. |^(a). Below this same temperature, it 
saturates at 0.1. Hence, the effective mass of the quasi- 
particle becomes ten times heavier than the bare elec- 
tron mass at low enough temperatures. The quasiparticle 
residue is almost isotropic throughout the Fermi surface. 
In Fig. ||(b) we present the drastic change of the quasipar- 
ticle dispersion (small filled circles) for T ~ 1/256 near 
the Fermi energy. The dashed curve is the noninteracting 
dispersion in the upper branch, (k) in Eq. ^, where the 
Fermi energy stays. The interacting quasiparticle disper- 
sion becomes dramatically flattened along the directions 
(0, 0) — (tt, tt) and (0, tt) — (0, 0), leading to a heavy effec- 
tive mass as well as large density of states at the Fermi 
energy. It is of interest to notice that the interacting en- 
ergy band appears separated into f electron (lower band) 
and c electron (upper band) dominating parts, which is 
absent in the noninteracting band. 

In order to understand the role of U in the one and two 
particle properties, we also performed the calculations for 
{/ = 3 and 4 with all the other parameters unchanged. In 
both cases, we found the indication of a ferromagnetic in- 
stability with decreasing temperature, which is signaled 
by a divergent behavior of XspiQ, 0) at q = (0, 0) in Fig. ||. 
Note that there is no finite temperature phase transition 
in two dimensions due to the Mermin- Wagner theorem. 
The inset clearly shows Xsp(0, 0) ~ exp{constant/T) at 
low temperatures. This magnetic instability happens be- 
cause the RKKY interaction between f electrons becomes 
dominating over the Kondo screening for large U ||l^. 
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FIG. 3. Static spin susceptibilities at g = (0, 0) for U = 3 
(filled circles) and ?7 = 4 (open circles). The inset shows the 
temperature dependence of Xsp(0, 0) in the low temperature 
regime. Notice the logarithmic scale for Xsp(0,0). 



In Fig. ^we present the density of states (Fig. ^(a)) and 
imaginary part of the self-energy (Fig. ^(b)) for [/ = 3, 
T = 1/2048 (solid curves) and C/ = 4, T = 1/290 (dot- 
ted curves). These temperatures are the lowest ones for 
each U where we found the convergent solutions for U^p 
in our numerical calculations. The inset in (a) shows the 
density of states in the wide range of frequency axis. Due 
to 2D strong spin fluctuations at low temperatures, the 
spectral weight is significantly suppressed in the low fre- 
quency regime (inset in Fig. ^(a)), compared with that 
for [/ = 2 (inset in Fig. |l](a)). To our surprise, a rela- 
tively sharp peak is still persisting at the Fermi energy 
even in the presence of strong spin fluctuations. In order 
to understand the nature of this peak, we computed the 
scattering rates in Fig. |4|(b). IniI]{kF,uj) shows a small 
dip at the Fermi energy in a large scattering background, 
which is different from the case for U — 2. This local min- 
imum of the scattering rates at the Fermi energy is re- 
sponsible for the sharp structure in the density of states. 
This feature comes from a partial screening of f electrons 
by conduction electrons for the following reason. In the 
critical region the imaginary part of the spin susceptibil- 
ity becomes singular like Irnxspi'l,'^) ^ S^i^ — fl) near 
q = (0, 0) for a ferromagnetic case, fl is the characteris- 
tic spin fluctuation energy which is less than 10~^ for the 
parameters used for U — 3 and 4. Hence, the imaginary 
part of the self-energy near the Fermi energy is given by 

ImT,f{kF,uj) J dv ImXsp{q,v)[B{v) + F{i> - w)] 

X I'mG^f{kp — q,ijj — v) 
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FIG. 4. (a) Density of states and (b) imaginary part of the 
self-energy at the noninteracting Fermi surface for e/ = 0.45 
and \/ = 1 with [/ = 3, T = 1/2048 (solid curves) and U = 4, 
T = 1/290 (dotted curves). The inset in (a) shows the density 
of states in the wide range of frequency axis. 
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FIG. 5. (a) Density of states and (b) imaginary part of the 
self-energy at the noninteracting Fermi surface for £/ = 0.45, 
;7 = 4 and T = 1/290 with increasing V. V = 1, 1.25 and 1.5 
are denoted as the solid, dotted and dashed curves, respec- 
tively. 
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Since ImXsp{<l,v) shows a nearly divergent behavior in 
the small region with the radius of approximately 7r/16 
from the origin in our calculations, the summation of q 
is over this area. Due to the mixture of f and c elec- 
trons through V , generally the hybridized band E^{k) 
has some dispersion or slope at the noninteracting Fermi 
surface. Because of this feature, the contributions from 
q ^ (0,0) give rise to larger scattering rates at 7^ 
than those at a; = 0, leading to a small dip at the Fermi 
energy in the scattering rates. This is in turn responsible 
for a relatively sharp peak in the density of states. Hence, 
a peak at w = in the density of states for [/ = 3 and 4 
is caused by mixing of f and c electrons, that is, a partial 
screening of local moments by conduction electrons. To 
confirm this argument, we also calculated the density of 
states and imaginary part of the self-energy by increasing 
the hybridization strength V with all the other parame- 
ters including the total electron concentration fixed. As 
mixing [V) is increased, the scattering rates from spin 
fluctuations are more suppressed near the Fermi energy 
in Fig. ^(b). The distance between the local maximum 
and the dip is given by the value E^^kp — q) at q = tt/16 
for all three different V. In a hypothetical situation where 
E'^{k) = constant (no mixing) and the same ImXspiQ, t^) 
is taken, the scattering rates would be maximum at the 
Fermi energy, leading to a pseudogap instead of a peak in 
the density of states. 

In summary, the formation of a coherent heavy-fermion 
state in the 2D periodic Anerson model has been stud- 
ied on the basis of the recently developed theory for the 



Hubbard model. For U — 2, the growing density of states 
and rapidly vanishing scattering rates near the Fermi en- 
ergy as well as linearly vanishing (in temperature) effec- 
tive magnetic moment below a characteristic tempera- 
ture, strongly support the Kondo screening of f electrons, 
leading to a heavy-fermion state. For J7 = 3 and 4, the 
dominance of RKKY interaction over Kondo screening 
at low temperatures indicates a magnetic instability at 
zero temperature. The persistence of a peak at cj = 
even in the presence of strong spin fluctuations is due to 
a partial screening of f electrons by conduction electrons. 
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